Vortices and Edge Reconstruction in Small Quantal Systems at High Angular 

Momenta 
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Vortices can form when finite quantal systems are set to rotate. In the limit of small particle 
numbers the vortex formation in a harmonically trapped fermion system, with repulsively 
interacting particles, shows similarities to the corresponding boson system, with vortices entering 
the rotating cloud for increasing rotation. We show that for a larger number of fermions, > 15, 
the fermion vortices compete and co-exist with (Chamon-Wen) edge-reconstructed ground states, 
forcing some ground states, for instance the central single vortex, into the spectrum of excited 
states. Experimentally, the fermion system could for instance be a semiconductor heterostructure, 
a quantum dot, and the corresponding boson system a magneto optical trap (MOT). 
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I. INTRODUCTION 

Macroscopically, the appearance of vortices is an every- 
day phenomenon. Hardly anyone has missed the swirl 
of water going down the drain, or the vortices forming 
while stirring in a cup of coffee. Correspondingly, on the 
microscopic level the formation of vortices and vortex lat- 
tices is known to lower the energy of a rotating quantum 
system. A vortex is then characterised by a singularity 
encircled by annular currents. 

In the 50's Abrikosov predicted the existence of vortices 
in superconductors subject to a magnetic field-. Al- 
most 50 years later, in 2001, Abo-Shaer et al? found 
vortices forming in rotating Bose-Einstein condensates. 
These many-body systems consist of paired electrons 
and bosonic atoms, respectively. Mean-field approaches 
are often employed to describe the complex structures 
that these systems comprise. For Bose-Einstein conden- 
sates, one often applies the Gross-Pitaevskii equation. In 
this way. Butts and Rokhsar'^ and later on Kavoulakis, 
Pethick and Mottelson"*'^ found successive transitions be- 
tween stable patterns of singly-quantised vortices, as the 
angular momentum was increased. 

Surprisingly vortex solutions, as discussed above, can also 
be found in small fermion systems with non-paired, repul- 
sively interacting particles at large angular momentaS^. 
Such systems could, for instance, be realized in quan- 
tum dota^ (small semiconductor islands confining a finite 
number of electrons) at high magnetic fields. Incredibly 
enough, the vortex formation by repulsively interacting 
fermions at high angular momenta is very similar to the 
bosonic case^: Irrespective of the system being bosonic 
or fermion, vortex solutions govern the ground states in 
the limit of large angular momenta. 
Setting a cloud of fermions rotating is equivalent to ap- 
plying a strong magnetic field. For an electron in a 
quantum dot, at a certain field strength the system spin- 
polarises, and forms a compact droplet where the N par- 
ticles occupy neighbouring single-particle orbitals with 
angular momenta m = 0, 1, 2, ....{N — 1), all belonging to 



the so-called lowest Landau level (LLL) with radial quan- 
tum numbers Ur — 0. This compact, so-called "maxi- 
mum density droplet" (MOD}- is the finite-size analog 
to integer filling factors in the bulk^". At larger fields or 
similarly, for higher total angular momenta, the droplet 
reconstructs its charge distribution. Different scenarios 
have been suggested (see& for a review): While in the 
low-iV limit, the formation of holes at the centre was dis- 
cussed by MacDonald, Yang and Johnson^'-'^^, Chamon 
and WerJ^ suggested that for larger droplets, the energy 
could be lowered by splitting off a (homogeneous) ring of 
electrons. In fact, it turned out later that the formation 
of a charge density wave along the edge was even more 
favourable^"^. 

In this article we re-investigate the reconstruction of the 
maximum density droplet (MDD) in the light of vortex 
patterns and the apparent analogy to the bosonic case. 
We show that, in agreement with the earlier findings of 
Yang et ali^, reconstruction of the MDD begins from 
the edge rather than the dot centre if the electron num- 
ber exceeds iV « 15. At smaller dot sizes, however, vor- 
tices enter the MDD from the edge, in much analogy to 
the bosonic case-. The universality between fermion and 
boson systems discussed by Toreblad et ali holds up to 
N « 15. For higher TV, not all vortex solutions appear 
as ground states. Instead edge reconstruction takes over 
at the onset of the MDD break-up. In the limit of very 
large angular momenta, vortices can co-exist with the 
usual (Chamon-Wen) edge-reconstructed ground states. 



II. MANY-BODY ASPECTS OF THE 
HARMONIC OSCILLATOR 



The Hamiltonian Hq, which we investigate, is of two- 
dimensional nature and describes an interacting finite 
many-body system under rotation (or, equivalently, sub- 
ject to a magnetic field, see appendix) which is confined 



2 



by a harmonic oscillator potential: 

i=l ^ ' i<j 

= H - QL, . (1) 

Here m is the particle mass and luq the confining fre- 
quency of the oscillator. This Hamiltonian has, in fact, 
become the standard choice to model finite quantal sys- 
tems - it equally well describes the physics of electrons 
trapped in a quantum do^, as bosonic or fermion atoms 
in a magneto-optical trap^. The two-dimensionality 
originates in the quantum dot from the semiconductor 
heterostructure. In the atom trap, it can be designed 
by the experimental set-up. In addition, strong rota- 
tional motion lowers the effective confinement perpen- 
dicular to the axis of rotation. This leads to a ground 
state with only the lowest energy level occupied along 
the z-directioniSii&. 

The operator —flLz adds rotation at an angular fre- 
quency Q. Since the total angular momentum, L, is a 
good quantum number to H^, it effectively subtracts a 
term JIL from the energy of the Hamiltonian H. As a 
result a larger will create a ground state with higher 
total angular momentum. Instead of seeking solutions at 
varying we follow the tradition in nuclear physics and 
solve H for a given L to obtain the total energy, , and 
the many-body state |^'_f,_b(L)). Note, that the lowest 
state for a fixed L is not necessarily the ground state of 
Hq. To find that state, we still need to subtract flL. 
As is usually done for these systems, the interaction be- 
tween the electrons in a quantum dot is modelled by the 
ordinary Coulomb repulsion: 

Viji - Tj) = ■ . 

A^Teeo\r^-rj\ 

For trapped, dilute atomic gases it is well known^"' that 
the interaction between the particles are more realisti- 
cally described by a zero-range potential of the form 
ad{ri — Tj), where a is the scattering lengtbi^. However, 
we have recently shown that the dominating features in 
the boson many-body spectra are to a large extent inde- 
pendent of the two-body interaction^, and very similar 
results were obtained both for short-range and Coulomb 
two-body forces in the lowest Landau level. For ease 
of comparison with the electron system under consider- 
ation, we thus restrict this study to the usual Coulomb 
interaction, both in the fermion and boson case. To sim- 
plify the comparison further, only one spin state will be 
considered: Spin-polarised fermions and bosons with zero 
spin. 

In second quantisation, here for fermions, we write 

OO oo 

I^(«l^li)clcj + 2 {^j\vcoui\kl)clc]ciCk, (2) 

where h stands for the single-particle part of the Hamilto- 
nian in Eq. {| is the harmonic oscillator basis 



and c] and Cj are the corresponding creation and anni- 
hilation operators. One can now simply construct the 
full Hamiltonian matrix and diagonalise it to obtain the 
energies and corresponding many-particle states. The 
many-particle solution for fermions is then given 

by an expansion in Slater determinants li^i), 

|vI/^)=^c,|V',). (3) 

i 

For the bosonic, symmetric solution \'^b), Eq. [21is simi- 
lar, with the fermion operators exchanged for the corre- 
sponding boson operators and 6^. 
This so-called "exact" diagonalisation is the brute-force 
method of many-body quantum mechanics. However, to 
name the method "exact" is misleading. A full Hamilto- 
nian matrix would be infinitely large, which of course 
would be impossible to handle. A further complica- 
tion is the fact that the size of the Hamiltonian ma- 
trix rapidly increases, making calculations for more than 
about six particles extremely time-consuming. In view 
of the structure of the single particle levels for high mag- 
netic fields^, we thus restrict the single particle states to 
the above mentioned lowest Landau level (LLL). Physi- 
cally this means that we ignore high energy excitations 
which would mix the lowest and higher Landau levels. 
For both bosons and fermions this is a common restric- 
tion, which is justifiable for dilute quantum gases at very 
high angular momenta, or at very high magnetic fields. 
It allows us to extend the studies to systems with up to 
25 particles. 



III. APPEARANCE OF VORTICES IN 
BOSONIC AND FERMIONIC SYSTEMS 

In weakly interacting bosonic systems, which are set ro- 
tating, successive transitions between stable vortex states 
have been found^ for increasing angular momentum. In 
particular, in a rotating and weakly interacting bosonic 
system with a short-range interaction, solutions of the 
Gross-Pitaevskii equation have revealed"*'^ a central one- 
vortex solution for the ratio of angular momentum to 
particle number being equal to one, / = Lfy/N = 1, a 
two-vortex solution at / = L^y/7V « 1.75 and a three- 
vortex state at I — Lfy/N 2.1 (in atomic units). 
Assuming similar vortex formation in the (unpaired) 
fermion system the vortex angular momenta are shifted 
by N{N — l)/2 as compared to the boson system. This 
shift corresponds to the many-body configuration of the 
MDD, which is the fermion equivalent to the boson many- 
particle ground state. An intuitive explanation is given 
by the configurations of the Fock states in the LLL, i.e. 
I 'n,onin2 ■ ■ ■ nm ■ ■ ■) , where is the number of par- 
ticles in the single-particle state with angular momen- 
tum TO. In the non-rotating ground state in the bosonic 
case all N spin-less particles are occupying the lowest 
single-particle orbital, the Fock state being of the form 
jA^OOOO • • • 000). For spin-polarised fermions, the Pauli 
principle demands single occupancies up to the Fermi 
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surface, and the Fock state of the MDD is a single Slater 
determinant of the form |111 • • • 1111000 • • -OOO). As the 
angular momentum of the MDD equals N{N — l)/2, a 
central single vortex in the fermion case should appear 
Lmdd)/N = 1, a two- vortex solution at 



at I — (^Ij-^Y 

I = {L^y'LMDD)lN 

correspondingly at I — 



1.75, and a three- vortex solution 
Lmdd)/N ^2.1. 



IV. COMPARISON BETWEEN ROTATING 
FERMION AND BOSON SYSTEM 

Let us now investigate the rotating fermion system as 
compared to the rotating boson system. In the lowest 
Landau level, at fixed angular momentum L the total 
energy is given by 



Er 



hujoiL + N) + (*F,s(i)|f |*F,s(i)), 



where the first term originates from the kinetic energy 
and the confining potential, and the second term from 
the interaction energy. An increase in hujQ{L + N) in- 
dicates a centre of mass excitation, whereas variations 
in {'i p,b{L)\v\'^ f,b{L)) imply changes in the interaction 
energy, which affects the configuration. 
Computationally, there now exist two more simplify- 
ing circumstances: (1) The kinetic energy and the con- 
fining potential constitute a diagonal term to the to- 
tal energy. It is therefore sufficient to diagonalise only 
{'^' F^BiL)\v\'^' f^b{L)) and add the term hujo{L + N) at 
the end to obtain the total energy. The outcome will be 
the same as including the diagonal terms in the diagonal- 
isation. (2) The interaction is independent of confining 
frequency, luq, since 



(^F,BiL)\v\^F,B{L))l 



V^{^f.b{L)\v\^f,b{L))\ 



This can easily be verified by writing the expression 
with second quantisation and comparing the integrals for 
= uj and luq = 1 a.u. The results presented are thus 
independent of confining frequency and it is only neces- 
sary to study loq — I a.u. as long as only the lowest 
Landau level is considered. 

The interaction energy versus angular momentum is dis- 
played in Fig. n for fcrmions (a) and bosons (b). Follow- 
ing the tradition in nuclear physics, the dashed line con- 
necting the lowest states at fixed L is called the "yrast" 
line — , and the spectrum of H the yrast spectrum. The 
interaction energy is plotted as a function of the above 
defined unit less variable /, which in a boson system 
corresponds to I — L/N and in a fermion system to 
I — [L ~ Lmdd)/N. By this construction we simphfy 
the comparison between boson and fermion systems, and 
it also becomes easier to compare systems of differing 
particle number. Note that Z = in the fermion sys- 
tem corresponds to the MDD. In the boson system, it 
describes the L — Q ground state. 

Not surprisingly, both in the boson and fermion spec- 
tra (as shown for = 12 in Figure the interaction 
energy decreases with increasing A larger Z-value im- 
plies a higher rotation, which pushes the particles further 
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. FIG. 1; Interaction energy versus angular momentum for a 
system of 12 harmonically trapped, spin-less particles inter- 
acting with a Coulomb repulsion, in (a) for fermions, and in 
(b) for bosons. The dashed line is the so-called yrast line, 
connecting the lowest-energy states. 



apart. However, the spectra reveal other striking similar- 
ities. Recalling the fact that, due to the Pauli principle, 
the fermion total angular momenta are displaced with 
respect to the boson systems by the angular momentum 
of the MDD, we realize that the spectra for fermions and 
bosons are very similar. The yrast line in both cases 
shows a pronounced kink at Z = 1. and fainter kinks are 
observed at I — 1.67 and at I — 1.83 (not shown in Fig.^). 
For the boson system we know these particular /-values to 
be related to the formation of vortices. According to the 
previously discussed results of the Gross-Pitaevskii equa- 
tions a one- vortex solution corresponds to exactly / = 1, 
and a two- vortex solution to / w 1.75. What happens in 
the fermion system? 

Before studying the many-particle states in more detail, 
it is interesting to explore which /-values correspond to 
ground states in the rotating system. In Figure El the 
ground states for increasing il have been found by sub- 
tracting D,L from El. Intriguingly, most of the /-values 
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FIG. 2: By taking El — for increasing total angular 



momentum the ground state values of I 



are found in a rotating system of 12 fermions as a function 
of the rotational frequency, Q.. For increasing rotation the 
ground state is shifted from the MDD {I = 0) to higher total 
angular momentum. In the figure successive transitions to 
states Z ~ 1, i ~ 1.6 and to I — 2.25 occur. For computational 
reasons higher Z-values have not considered. 



corresponding to kinks of the yrast line, including inter- 
mediate values I — 0.92 and 1.5, now reappear. (Note 
that to increase above 1 a.u. has no physical mean- 
ing. If the rotational frequency exceeds the confining 
frequency the effective confinement disappears, see ap- 
pendix.) 

To extract information from the many-particle wavefunc- 
tions, one must process the data in some way. Two in- 
formative quantities are the occupancies, Pm, of the sin- 
gle particle states, which illustrates to what extent each 
single particle state (in the LLL) is present in the many- 
body configuration, and the radial particle densities, 



JV 



F,B 



where r = (r^O) in polar coordinates. The occupancies 
describe which single particle states dominate. If, for in- 
stance, Pm=Q has a magnitude close to zero the radial 
density will be very small at r = 0, since the single par- 
ticle state (in the LLL) with m — has its maximum at 
r = 0. Similarly a minimum at Pm=i implies a slightly 
off-centred minimum in the radial density. Hence, a cen- 
tral single vortex should appear as a minimum for m = 
in the occupancy and at r = in the radial density, and 
a two-vortex state as an off-centred minimum both in 
the occupancy and in the radial density. In Figure 13 
the occupancies are plotted for the ground state val- 
ues I — 0.92, 1, 1.5, 1.67 and 2.25, and the excited state 
I = 1.83. The insets show the radial densities. The last 
state has been included for completeness. If we compare 
with the ground states in the bosonic system, these are 






5 10 15 20 
angular momentum, m 

FIG. 3: The occupancies of the single particle states for a 
system of iV = 12 fermions. The values / = 0.92, 1, 1.5, 1.67 
and 2.25 correspond to the ground state values in a system 
rotating at an increasing frequency. The excited state I — 
1.83 has been included for completeness. The insets show the 
radial densities from r — —Ga*g to r — 6a*g for each /-value. 



I = 1, 1.67, 1.83, and 2.25. Even if the state I = 1.83 in 
the fcrmion system is not a ground state in the rotat- 
ing system, it is the lowest-energy state at that particu- 
lar I in the yrast spectrum. In the boson case we know 
that a central single vortex appears at / = 1 and a two- 
vortex solution is found at Z « 1.75. As expected, in the 
fermion system, the central single vortex is clearly visible 
at I = 1. The two-vortex state, however, is more diffi- 
cult to identify. In contrast to mean field methods, for 
an exact calculation with good angular momentum the 
particle density has circular symmetry and thus does not 
display the internal structure directly. The vortices are 
averaged over the entire circle and only show as two dips 
in the radial densities. The ground state with I = 1.67 
and the excited state with I = 1.83 showing two minima 
in the radial density thus constitute likely candidates for 
a two-vortex state. 

Another implication of many-vortex solutions in the ex- 
act result can be found in the occupancies, which display 
some interesting systematics. To understand these sys- 
tematics better we chose to study only the structure of 
the Fock state with the largest magnitude in the expan- 
sion of l^piL)) (see Eq.^. 

In Table^the dominating configurations and their ampli- 
tudes are shown for ground state Z-values and the excited 
state I = 1.83 in the case of iV = 12. For I = 0.92 a single 
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0.44 
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1.83* 


0.55 


1001111111111100 


2.25 


0.28 


1110011111111100 



TABLE I: The dominating many-particle configurations for 
A'^ = 12 for tiie ground state i-values. Note the single, re- 
spectively double, hole for I — 1 and I — 1.83, which are the 
expected central single and double vortex states. {I — 1.83 is 
an excited state, as marked by a star). 



hole has entered the MDD. This corresponds to the mini- 
mum shown in the occupancy in Figure|21 At ^ = 1 it has 
reached the centre of the cloud and the single hole state 
has become very dominant in the many-body configura- 
tion (consider its magnitude). Correspondingly, the ra- 
dial density shows a clear central single vortex. At ^ = 1.5 
a double hole starts to enter the MDD. It moves further 
towards the centre at I — 1.67 and in the excited state 
I = 1.83 it has reached the position, where its amplitude 
is the largest. This configuration seems to be a double 
hole close to the centre of the fermion cloud. Could it be 
a two- vortex state? At I = 2.25 a triple hole begins to 
enter the MDD. 

As compared to mean-field approaches, these "exact" 
solutions retain the circular symmetry of the Hamilto- 
nian. In order to reveal the internal structure of the 
many-body wave function, methods such as calculating 
the pair correlation, looking at a slightly asymmetric 
confinement potentiali^ or adding a small perturbation? 
have to be used. Especially for small N, the above 
methods have shown useful: When, for example, study- 
ing pair correlations, fixing one of the very few particles 
was found to cause too much disturbance in the fermion 
many-particle configuration to get a meaningful result. 



V. PAIR CORRELATIONS 

The pair correlation fixes one particle at a point Ta and 
calculates the density of the remaining iV — 1 particles. 

For large bosonic systems this shows the vortices clearly, 
see for instance Kavoulakis et a&. 

However, in the limit of small N, particularly in the 
fermion case, the removal of one particle effects the den- 
sity of the remaining particles significantly, and the re- 
sult can be ambiguous: Due to the Pauli principle the 
point, ta, at which one fermion has been fixed, will be 
avoided by the other fermions. A minimum will arise in 
the pair correlation at this point called the exchange hole. 



This makes the detection of vortices in the pair correla- 
tions very tedious. Particularly in the limit of small N, 
a small displacement of the reference point alters the 
pair correlations considerably. 

To come about this difficulty, two schemes have pre- 
viously been suggested: (i) adding a perturbation to 
the harmonic potential, and (ii) breaking the rotational 
symmetry, for instance by deforming the confinemen^iS. 
However, it appears that TV = 12 is sufficiently large 
to study the pair correlations. The pair correlations for 
= 12 have been plotted in Figure 01 (fermions) and 
Figure |S1 (bosons) for the ground state values of I in both 
the fermion and bosonic case. For completeness, the ex- 
cited state I = 1.83 has also been added in the fermion 
case. Apart from this state and the two intermediate 
states I = 0.92 and 1.5, there is an exact correspondence 
in the ground state ^-values of the fermion and boson 
systems. With respect to angular momentum, this cor- 
responds to the shift of the angular momentum Lmdd 
of the maximum density droplet in the fermion system. 
Studying the pair correlations for these ^-values, the en- 
tering vortex in the rotating system can clearly be seen, 
first changing to a central single vortex state at Z = 1, 
and then to a double at I = 1.83. This (excited) state 
is the same two-vortex solution as in the bosonic case, 
appearing at similar I. The solution at / — 2.25 shows a 
third vortex entering, as expected from the occupancies. 
Note how the Pauli principle affects the fermion pair cor- 
relations by introducing the exchange hole (as marked by 
a black dot in Fig. 0)). Naturally, this exchange hole is 
not present in the corresponding boson pair correlations. 
These are plotted in Figure [S] for comparison. We note 
that although the pair correlations in general are a bit 
ambiguous and the results depend strongly on the choice 
of the reference point, this method was earlier shown^ to 
yield a similar picture of vortex entry into the rotating 
cloud as the Gross-Pitaevskii approacl^. 
Compared to the mean-field picture the minima will not 
equal zero due to the zero-point oscillations present in 
the exact diagonalisation results. 



VI. EDGE RECONSTRUCTION 

For particle numbers larger than approximately 15 the 
central single vortex disappears from the lowest state at 
I = 1 and moves up into the spectrum of excited states. 
It appears that, in the fermion system, the vortex con- 
figuration is not as favourable anymore-^ 
In figure El the yrast spectrum for = 20 fermions is 
displayed for increasing values of /. A second order poly- 
nomial, aP + bl + c, has been subtracted from the in- 
teraction energy in order to display the structure of the 
spectrum in greater detail. In fact, below I = 0.9, the 
entering vortex state dominates as it did for smaller par- 
ticle numbers. However, at I — 0.9, there exists a level 
crossing in the yrast spectrum. The state with the enter- 
ing vortex becomes an excited state, which corresponds 
to the central one- vortex configuration at Z = 1, while a 
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FIG. 4: Contours of the pair correlation (left panel; the black 
dots mark the exchange hole and the scale goes from minimum 
at dark blue to maximum at dark red) and radial density pro- 
files (right panel) for iV = 12 at Z = 0.92, 1, 1.5, 1.67, 1.83 and 
2.25. The vortices seem to enter the cloud from the outer 
regions, and then move towards the centre. At Z = 1 a central 
single vortex can be seen and at I — 1.83 a two- vortex solution 
exists. 



state with a distinct edge becomes most prominent above 
I = 0.9. Similarly, the Z-values at which the actual ground 
states occur in the rotating spectrum, are altered. The 
central one-vortex, respectively the entering two-vortex 
at Z = 1 and / « 1.7, which we saw for 12 fermions, no 
longer appear as ground states. Instead, configurations 
around I — 0.6, I — 1.2 and I — 1.8 have become more 
favourable in energy. All three constitute states with a 
clear edge, and in particular / = 0.6 and I = 1.2 are in- 
teresting since their edges consist of localised fermions. 




FIG. 5: Contours of the pair correlation for A'^ = 12 bosons, 
plotted as in Fig. 2] The reference point is marked with a 
black dot. 



so-called Chamon-Wen edges (see Figure 01. This is in 
contrast to the rotating boson system which is known^ 
to have a central one- vortex ground state at Z = 1. 



FERMIONS 




FIG. 6: Low-lying many-body states as a function of / for N = 
20 fermions. A polynomial al^ + hi + c has been subtracted 
from the exact many-body energies to display the structure 
of the spectrum in greater detail. Just below the central unit 
vortex, a level crossing occurs: The single vortex state (green 
line) moves upwards in energy, and is replaced by an edge 
reconstructed ground state. 

A Chamon-Wen edge becomes favourable when a higher 
rotation (or magnetic field) effectively decreases the ra- 
dius of the single particle states as compared to the con- 
fining trap. For a rotating system at high angular mo- 
mentum this would imply an increase in density rela- 
tive to the confinement with increasing rotational mo- 
tion, hence effectively compressing the system. The re- 
pulsive Coulomb interaction strongly inhibits this, and 
the maximum density droplet in effect redistributes its 
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Rotational frequency, Q [a.u.] 



FIG. 7; Ground state i-values for increasing rotation (fi) no 
longer contains the central single and double vortex states in 
the 20 fermion system. Instead edge reconstructed states, like 
I = 0.6 and I — 1.2 become dominating. 



particles from lower to higher single particle momentum, 
leaving some states with m < N —1 unoccupied, giving a 
ground state with L > N{N —l)/^^^. Figure |H1 shows the 
occupancies and radial densities of some of the lowest en- 
ergy states along the yrast line for increasing I. Both the 
occupancies and the radial densities imply a Chamon- 
Wen edge, i.e. that there exists an edge of redistributed 
fermions and an inner dense MDD droplet. Also visible is 
the entering single vortex configuration at I = 0.9 which 
becomes the second excited state in the yrast spectrum 
at ^ = 1 (i = 3). This is the central one- vortex state. 
Studying also the pair correlations for = 20 pictured 
in Figure El for I = 0.6 respectively I = 1.2, the edge 
of localised fermions is clearly visible. These pair cor- 
relations compare well with the results by Reimann et 
al^, showing electron localisation along the Chamon- 
Wen edge within mean field current spin density func- 
tional theory (CSDFT). Note that in the exact result 
the localisation becomes less pronounced due to the zero 
point oscillations. At the first glance one might draw the 
conclusion that the central single vortex disappears from 
the lowest energy at / = 1 since we have increased the 
number of fermions in the trap. This suggests a higher 
particle density, which does not allow for a central mini- 
mum, the single vortex. However, the interaction energy 
only scales with a factor depending on ujq (see Sec. II). 
Hence, our results show that the central single vortex 
state will no longer be favourable irrespective of confine- 
ment frequency. 

Solutions where vortices enter at the edge of the droplet, 
can co-exist with edge reconstructed configurations. The 
entering second vortex, as well as the (now off'-centred) 
first vortex appear as two clear minima. Consider the 
pair correlation for I — 1.2 in Figure |51 A more pro- 
nounced entering two- vortex state shows in the pair cor- 
relation of I — 1.5 along the yrast line, also in Figure |51 
This agrees well with the double hole in the occupancy 
of Figure |S1 for the same state. The entering second vor- 
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FIG. 8: Occupancies of the single particle states for a system 
of A'' = 20 fermions, shown for the lowest states for increasing 
I. Corresponding radial densities from r = —7a*g to r — 7ag 
are shown in the insets. The central single vortex state has 
now disappeared from the lowest state (i = 1) and has become 
the 2nd excited state (i = 3) at / = 1. 



tex can also be seen in the energy spectrum in Fig. |^ 
as oscillations with a period of A/ — 0.1. As explained 
earlier, by subtracting a second order polynomial, the 
oscillations in the yrast spectrum are enhanced. In Man- 
nmenet ali^ this was discussed in greater detail. It was 
found that the oscillations change period at higher I, first 
to become three-fold and then to become four-fold. In 
Figure ITUl we compare the oscillation periods with holes 
in the Fock states, by studying the dominating configu- 
rations for increasing I. For simplicity, the lowest energy 
at each /-value {EL(i = 1)) have been subtracted, i.e. 
we study the excitation spectrum. The (red) circles de- 
scribe the entering vortex state, |11 • • • 11011 • • • 10 • • •), 
which is seen to become the second excited state in the 
yrast spectrum at / = 1. Earlier, this was found to be 
the central single vortex configuration. The (green) tri- 
angles are then seen to move downwards in the spectrum 
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FIG. 9: Pair correlations for N=20 fermions for some values 
of I along the yrast line. In addition, I = 0.6 and I = 1.2 
are ground states in the rotating system. Especially note 
the clear Chamon-Wen edge of the I — 0.6 state, which can 
be compared to the CSDFT result in Reimann et al^ with 
I = 0.75. 



spectrum as oscillations with period Al = 0.1. Moving 
the double hole (00) one step to the left implies an in- 
crease in total angular momentum of two. The entering 
double vortex, thus, only appears at every second L (in 
the Figure Al = 0.1), corresponding to the double os- 
cillations. Similarly, the entering triple vortex, (black) 
diamonds, with Fock state |11 • • • 1100011 • • • 10 • • •), can 
only occur at every third L, corresponding to a three- 
fold oscillation period. This is in exact agreement with 
Manninen et ali^. 



VII. SUMMARY AND CONCLUSIONS 

To conclude, the results presented have shown that there 
indeed exists some universal features between rotating 
systems of bosons and (unpaired) fermions. All calcula- 
tions have been performed using the method of exact di- 
agonalisation, and the wavefunction data have been eval- 
uated with the aid of occupancies and pair correlations. 
Especially, the tendency of these systems to form vortices 
at similar angular momenta have been explored. How- 
ever, for larger particle numbers than approximately 15 
the central single vortex state becomes an excited state 
and the ground state, instead, lowers its energy by re- 
constructing to form an edge of localised fermions. This 
is in agreement with earlier mean-field computations, see 
for instance Reimann et eli^. 
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APPENDIX: ANALOGY BETWEEN ROTA- 
TION AND MAGNETIC FIELD 

Consider the Hamiltonian of the rotating system, Eq. 
in Sect. II. Re-writing Hq, we show that it equally well 
corresponds to a harmonic oscillator subject to an effec- 
tive magnetic field of strength B — W x {milz x f/e) = 
(2rnf2/e)z and with a confinement frequency (wq — r2^)i&: 



FIG. 10: Yrast spectrum for A'' = 20, with the lowest energy 
at fixed / {EL{i ~ 1)) subtracted. The (red) circles mark a sin- 
gle zero entering the Fock state of the MDD, the (green) trian- 
gles a double hole, for example | 11 ■ ■ • 11100111 • • • 11000 ■ ■ ■) 
and (black) diamonds a triple hole. 
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for increasing I. These are the entering double vortex 
states, |11 • • • 110011 • • • 10 • • •), appearing in the energy 



As a consequence, the rotation effectively lowers the con- 
finement in x-y-direction to (wg — 51^). 
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